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ABSTRACT 



Context. Two classes of gamma-ray bursts have been identified in the BATSE catalogs characterized by durations 
shorter and longer than about 2 seconds. There are, however, some indications for the existence of a third class. Swift 
satellite detectors have different spectral sensitivity than pre-Swift ones for gamma-ray bursts. Therefore we reanalyze 
the durations and their distribution and also the classification of GRBs. 

Aims. We analyze the bursts duration distribution, published in The First BAT Catalog, whether it contains two, three 
or more groups. 

Methods. Using The First BAT Catalog the maximum likelihood estimation was used to analyze the duration distri- 
bution of GRBs. 

Results. The three log-normal fit is significantly (99.54% probability) better than the two for the duration distribution. 
Monte-Carlo simulations also confirm this probability (99.2%). Similarly, in previous results we found that the fourth 
component is not needed. The relative frequencies of the distribution of the groups are 7% short 35% intermediate and 
58% long. 

Conclusions. Similarly to the BATSE data, three components are needed to explain the BAT GRBs' duration distri- 
bution. Although the relative frequencies of the groups are different than in the BATSE GRB sample, the difference 
in the instrument spectral sensitivities can explain this bias. This means theoretical models may be needed to explain 
three different type of gamma-ray bursts. 
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1. Introduction 

It has be en a great challenge to cl assify gamma-ray b ursts 
(GRBs). iMazets et all (| 19811 ) and iNorris et alJ (| 19841 ) sug- 
gested there might be a separation in their duration dis- 
tribution. Using The First BATSE Catalog, Kouveliotou 
et al. (1993) found a bimodality in the distribution of the 
logarithms of the durations. In that paper they used the pa- 
rameter Tgp (the time in which 90 % of the fluence is accu- 
mulated (Kouveliotou ct al., 199 3)) to characterize the du- 
ratio n of G RBs dMcBreen et all. I1994 iK oshut et all 119961: 



In a previous paper using the Third BATSE Catalog 
(jMeegan et all . 119961) iHorvatbl (|1998l ) showed that the du- 
ration (T 90 ) distribution of GRBs observed by BATSE 
could be well fitted by a sum of three log-normal distri- 
butions. We find it statistically unlikely (with a probability 
~ 10~ 4 ) that there are only two groups. Simultaneously, 
Mukherjee et al. (1998) report the finding (in a multidimen- 
sional parameter space) of a very similar g roup structure 
of GRBs. Somewhat later, s everal authors (jH akkila et al 



2000; lBalastegui et all l200lURaianiemi fc Mahonenl. 12 002: 



iBellil . 119971 : iPendleton et all Il997l) . Today it is widely ac 



Hakkila ct al 



cepted that the physics of these two groups (also called 
"subclasses" or simply "classes") are differe nt, and these 
two kinds of GRBs are different phenomen a (|Norris et all 
[20011: iBalazs et all [20031 : iFox et all 12001 . In the Swift 
database the measured redshift distribution for the two 
grou ps are also different, for s hort bursts the median is 
0.4 (lO'Shaughnessv et all I2008D and for the long ones it is 
2.4 jBagolv et all l2006f ) . 

The bimoda l distribution was furth er quantified in 
another paper (jKouveliotou et all 119961 ). where a two- 
log-normal fit wa s made; the best par amet ers of the fit 
were p ublished in lMcBreen et al.l (|1994[ ) and lKoshut et al.l 
(fl99T 



2003; 



120041: Hakkil a et al 



Borgonovol t 

12004 : Chattopadhvav et all 20071) included more phys- 
ical parameters in the analysis of the bursts (e.g. 
peak-fluxes, fluences, hardness ratios, etc.). A clus- 
ter analysis in this multidimensional parameter space 
suggests the ex istence of the t hird (" intermediate") 



as well dMukheriee et aT . 19981: I Hakkila et al 
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group 

| 2000t [Balastegui et all 120011: iRaianiemi k Mahonenl . l2002t 
IChattopadhvav et al.l . I2007T ). The physical existence of the 
third group is, however, still not convincingly proven. 
However, the celestial distrib u tion of the third group 
is anisotropic flMeszaros et all 120001 : iLitvin et all l200lh 
Mag liocchetti et all 120031 ). All these results mean that the 
existence of the third intermediate group in the BATSE 
sample is acceptable, but its physical meaning, importance 
and origin is less clear than those of the other groups. 
Hence, it is worth studying new samples if their size is large 
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Fig. 1. Duration distribution of the observed BAT bursts. 



enough for statistical an alysis. In the HETE-II database 
(jVanderspek et al.l . 120041) the re are only 104 GRBs and in 
the Swift first BAT database (|Sakamoto et all 120081) there 
are 237 GRBs. Therefore, in this paper we use the Swift 
data because of its better statistics. 

In Sect. 2 we discuss the method used in the paper. In 
Sect. 3 uni-, bi-, tri- and tetra- modal log- normal fits made 
by using the maximum likelihood method are discussed. In 
Sect. 4 one thousand Monte-Carlo simulations are shown 
investigating the significance of the fits. In Sect. 5 we dis- 
cuss some further details. The conclusions are given in Sect. 
6. 



2. The method 

There are several methods to test significance. Fo r example 
the y 2 method which we used in our first paper (Horvath, 
119981) to analyze the T go distribution of the BATSE bursts 
is not useful here, because of the small population of short 
bursts in the Swift sample. 

In the Swift BAT Catalog ([Sakamoto et al.l . l2008f) there 
are 237 GRBs, of which 222 have duration information. Fig. 
1. shows the log Tgo distribution. To use the % 2 method one 
has to bin the data. If the number of counts within some 
bins is small the method is not applicable. The Maximum 
Likelihood (ML) method is not sensitive to this problem, 
therefore for the (Swift) BAT bursts the maximum likeli- 
hood method is much more appropriate. 

The ML method assumes that the probability density 
function of an x observable variable is given in the form of 
g(x,pi, ...,pk) where p\,...,pk are parameters of unknown 
value. Having N observations of x one can define the like- 
lihood function in the following form: 



N 



I = Y[g(xi,pi, -,Pk), 



or in logarithmic form (the logarithmic form is more con- 
venient for calculations): 



A' 



L = log/ = ^log(g(xi,pi, ...,p k )) 



Table 1. The best parameters for the two log- normal fit of 
the GRB duration distribution. 





Duration(logTgo) 


a(logT 90 ) w 


short 


-0.456 


0.501 16.3 


long 


1.606 


0.507 205.7 
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Fig. 2. Fit with two log-normal component for the duration 
distribution of BAT bursts. 



The ML procedure maximizes L according to pi, ...,pk- 
Since the logarithmic function is monotonic the logarithm 
reaches the maximum where I does as well. The confidence 
region of the estimated parameters is given by the following 
formula, where L max is the maximum value of the likelihood 
function and L is the likelihood function a t the true value 
of the parameters ([Kendall k, Stuartl . I1973T ): 



2(L r 



Lo) 



(3) 



3. Log-normal fits of the duration distribution 

Similarly to iHorvathl (|2002t) we fit the log Tgo distribution 
using ML with a superposition of k log-normal components, 
each of them having 2 unknown parameters to be fitted with 
N = 222 measured points in our case. Our goal is to find 
the minimum value of k suitable to fit the observed distri- 
bution. Assuming a weighted superposition of k log-normal 
distributions one has to maximize the following likelihood 
function: 




y^ y Wifl{xi,logTi,ai) 



(4) 



where wi is a weight, /; a log- normal function with logT; 
(1) mean and oi standard deviation having the form of 



QE-logT,)' 



and due to a normalization condition 

k 



(2) J2 W[=N - 



(5) 



(6) 



I. Horvath et al.: Classification of Swift's gamma- ray bursts 



3 



35 




-1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 
log T90 

Fig. 3. Fit with three log-normal component for the dura- 
tion distribution of BAT bursts. 



Table 2. The best parameters for the three log-normal fit 
of the GRB duration distribution. 





Duration(logTgo) 


a(logT 9Q ) 


w 


short 


-0.473 


0.48 


16.2 


long 


1.903 


0.32 


129.1 


intermediate 


1.107 


0.35 


76.7 



Table 3. The improvement of the likelihood and the sig- 
nificancies. 



i 


Limax 


Limax ^(i-l)max 


V 


2 


983.317 






3 


989.822 


6.505 


0.9954 


4 


990.323 


0.501 


0.200 



We used a simple C++ code to find the maximum of 
Lfc. Assuming only one log- normal component the fit gives 
Limax — 951.666 but in the case of k=2 one gets L<x max — 
983.317 with the parameters given in Table 1 and the solu- 
tion displayed in Fig. 2. 

Based on Eq. ([3]) we can infer whether the addition of a 
further log-normal component is necessary to signifincantly 
improve the fit. We take the null hypothesis that we have 
already reached the the true value of k. Adding a new com- 
ponent, i.e. moving from k to k + 1, the ML solution of 
Lkmax change to La.+x)maxi but Lq remained the same. In 
the meantime we increased the number of parameters by 
3 (wk+i, logT k+1 and <J( k +i))- Applying Eq. © to both 
Lkmax and L( k +i)max we get after subtraction 

2{L(k+l)max ~ Lkmax) ~ x\ ■ (?) 

For k = 1 Li max is greater than L\ max by more than 30, 
which gives for \\ an extremely low probability of 5.88 x 
10~ 13 . This means that the two log- normal fit is really a 
better approximation for the duration distribution of GRBs 
than one log-normal. 

Thirdly, a three-log-normal fit was made combining 
three fk functions with eight parameters (three means, 
three standard deviations and two weights). For the best 
fit parameters see Table 2. The highest value of the log- 
arithm of the likelihood (L^ max ) is 989.822. For two log- 
normal functions the maximum was Limax = 983.317. The 
maximum thus improved by 6.505. Twice this is 13.01 which 
gives us the probability of 0.461% for the difference between 
Limax and L^ max is being only by chance. Therefore there 
is only a small chance the third log-normal is not needed. 
Thus, the three-log-normal fit (see Figure 3.) is better and 
there is a 0.0046 probability that it was caused only by 
statistical fluctuation. 

One should also calculate the likelihood for four log- 
normal functions. The best logarithm of the ML is 990.323. 
It is bigger by 0.501 than it was with three log-normal func- 
tions. This gives us a low significance (80.1%), therefore the 
fourth component is not needed. In Table 3. we summarize 
the improvement of the likelihood and the probabilities they 
give us. 



4. 1000 Monte-Carlo simulations using the 
two-component fit 

We can check the 0.0046 probability, which we get for 
the maximum likelihood calculation, using a Monte-Carlo 
(MC) simulation and adopting the following procedure. 
Take the two-log-normal distribution with the best fitted 
parameters of the observed data, and generate 222 num- 
bers for Tgo whose distribution follows the two-log-normal 
distribution. Then find the best likelihood with five free 
parameters (two means, two dispersions and two weights; 
but the sum of the last two must be 222). Next we per- 
form a fit with the three-log-normal distribution (eight free 
parameters, three means, three dispersions and two inde- 
pendent weights). Finally, we take the difference between 
the two logarithms of the maximum likelihoods that gave 
one number in our MC simulation. 

We have carried out this procedure for 1000 simulations 
each with 222 simulated logTgoS. There were 8 cases when 
the log-likelihood difference was more than the one obtained 
for the BAT data (6.505). Therefore the MC simulations 
confirm the result obtained by applying Eq. |J7J| and give 
a similar (0.8%) probability that a third group is merely a 
statistical fluctuation. 



5. Discussion 

It is possible that the fit using three log-normal functions 
is accidental, and that there are only two types of GRBs. 
However, the probability that the third component is only 
a statistical fluctuation is 0.5-0.8 %. 

One can compare the burst group weights with pre- 
vious re sults. BAT sensitivity is different to BATSE sen- 
sitivity (|Fishman et all 11994 iBandl . l2003f) . BAT is more 
sensitive at low energies which means it can observe more 
X-ray flashes and soft bursts and probably fails to de- 
tect many hard bursts (typically short ones). Therefore 
one expects more long and intermediate bursts and fewer 
short GRBs. In the BAT data set there are only a few 
short bursts. Our analysis could only find 16 short bursts 
(7%). The robustness of the ML method is demonstrated 
here because a group with only 7% weight is identified. 
Previously in the BATSE database intermediate bursts 
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Table 4. Mean hardness ratios and their standard errors. 



type 



logH43 std.err.log H43 logH32 std.err.logH32 



short 
intermediate 
long 



-0.023 
-0.185 
-0.092 



0.023 
0.014 
0.011 



0.276 
0.112 
0.184 



0.030 
0.018 
0.012 



were identified by many research groups. However, in this 



of BATSE GRBs (Mukheriee et all Il998t lHakkila et al. 


2000; Balastegr 


d et all l200lHRaianiemi & Mahonenl. 2002 


Horvath ct al., 


200l. 



We calculated the mean logH43 and logH32 hardness 
ratios and the standard errors of the groups as given in 
Table 4. The table demonstrates, as expected from previ- 
ous studies, that short bursts are the hardest and the in- 
termediate duration group is the softest among the GRBs 
detected by Swift. Using a t-test ([Kendall fc StuartJ . Il973j ) 
the intermediate bursts differ very significantly (99.9% in 
both hardnesses) from the other two and the short and the 
long bursts are also significantly different in logH43 (98%) 
an d loqH32 (9 9 %). 

ICline et al l (|2005l ) claims that GRBs with T 90 < 0.1 s 
form a separate group. Out of the 222 Swift GRBs used 
in our analysis only 4 (25 % of the short population) have 
a duration of Tqq < 0.1 sec. This very low frequency of 
these very short GRBs does not allow a detailed statistical 
analysis. Additionally, the Swift satellite is less sensitive to 
these types of bursts due to trigger criteria (|McLean et all 
2004) 



6. Conclusions 

1. Assuming that the T90 distribution of the short and 
long GRBs is log-normal, the probability that the third 
group is a chance occurance is about 0.5-0.8 %. 

2. Although the statistics indicate that a third component 
is present, the physical existence of the third group is 
still debatable. The sky distribut ion of the third compo - 
nent is anisotropi c as pr oven by llvleszaros et al.l (fe OOO) 
and Litvin et al] (|200ll ). Alternatively lHakkila etlrt] 
( 2000) believe the third statistically proven subgroup is 
only a deviation caused by complicated instrumental ef- 
fects, which can reduce the duration of some faint long 
bursts. This paper does not deal with this particular 
effect, however the previously studied BATSE sample 
shows a similar group structure. This agreement sug- 
gests that the third component is possibly real, not an 
instrumental effect (the BATSE detectors and the Swift 
BAT arc different kinds of instruments). 

3. The observed frequencies in the three classes are differ- 
ent for BATSE and BAT. Both samples are dominated, 
however, by the long bursts. The short bursts are less 
populated in BAT than in BATSE but the intermediate 
group is more numerous. This is understandable, since 
BAT is less sensitive in high energy than BATSE was 
and more sensitive in low energy and short bursts are 
the hardest group and intermediate ones are the softest. 
Therefore BAT can observe more intermediate bursts 
and much fewer short ones than BATSE did. 



4. The existence and physical properties of the intermedi- 
ate group need further discussion to elucidate the reality 
and properties of this class of GRBs. 
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